The fc-th derivatives of the immanant and the 
X-symmetric power of an operator 

Sonia Carvalho* Pedro J. Freitas^ 
February, 2013 

Abstract 

In recent papers, R. Bhatia, T. Jain and P. Grover obtained for- 
mulas for directional derivatives, of all orders, of the determinant, the 
permanent, the m-th compound map and the m-th induced power 
map. In this paper we generalize these results for immanants and for 
other symmetric powers of a matrix. 

1 Introduction 

There is a formula for the derivative of the determinant map on the space of 
the square matrices of order n, known as the Jacobi formula, which has been 
well known for a long time. In recent work, T. Jain and R. Bhatia derived 
formulas for higher order derivatives of the determinant ([?]) and T. Jain 
also had derived formulas for all the orders of derivatives for the map A m 
that takes an n x n matrix to its m-th compound QlOj). Later, P. Grover, 
in the same spirit of Jain's work, did the same for the permanent map and 
the for the map V m that takes an n x n matrix to its m-th induced power. 
The mentioned authors extended the theory in [BJ. 
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This paper follows along the lines of this work. It is known that the deter- 
minant map and the permanent map are special cases of a more generalized 
map, which is the immanant, and the compound and the induced power of a 
matrix are also generalized by other symmetric powers, related to symmetric 
classes of tensors. These will be our objects of study. 

2 Immanant 

We will write M„(C) to represent the vector space of the square matrices 
of order n with complex entries. Let A e M n (C) and \ be an irreducible 
character of C. We define the immanant of A as: 

n 

creSn 1=1 

In other words, d x : M„(C) — > C is a map taking an n x n matrix to its 
immanant. This is a differentiable map. For X G M n (C), we denote by 
Dd x A(X) the directional derivative of d x at A in the direction of X. 

We denote by A(i\j) the n x n square matrix that is obtained from A by 
replacing the i-th row and j-th column with zero entries, except entry (i, j) 
which we set to 1. We define the immanantal adjoint adj x (v4) as the n x n 
matrix in which the entry (i, j) is d x (A(i\j)). This agrees with the definition 
of permanental adjoint in [131 P- 237], but not with the usual adjugate matrix 
(we would need to consider the transpose in that case). This is a matter of 
convention. 

We obtain the following result similar to the Jacobi formula for the de- 
terminant. 

Theorem 2.1. For each X e Af„(C), 

Dd x (A)(X) = tr(adj x (A) T X). 

Proof. For each 1 < j ' < n, let A(j;X) be the matrix obtained from A by 
replacing the j-th column of A by the j-th column of X and keeping the rest 
of the columns unchanged. Then the given equality can be restated as 

n 

Dd x (A)(x) = J2d x (Mr,x)). (1) 

i=i 
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On the other hand we note that Dd x (A)(X) is the coefficient of t in the 
polynomial d x (A + tX). Using the fact that the immanant is a multilinear 
function of the columns we obtain the desired result. □ 

Again using the fact that the immanant is a multilinear function, we 
notice that for any 1 < i < n, we have 

n 

rf x(^) = X^AWb'))- 

i=l 

Using this and (QQ), we get that 

n n 

Dd x (A)(X) =^^x ij( i x (i(i|j)). (2) 

i=l j=l 

We will generalize the expressions ([I]) and ([2]) for the derivatives of all 
orders of the immanant. 

We now turn to derivatives. Let Vi, . . . , V n be n vector spaces over C, 
and let (j) : V\ X . . . x V n — > C be a multilinear form. For A, X 1 , . . . X k G 
Vi, ... , V n , the k-th derivative of <fi at A in the direction of (X 1 , . . . , X k ) is 
given by the expression 

^A + hX 1 + ...+t k X k ). 

tl =...=t k =o 

This is a multilinear function of M k (C). In particular, if we consider 
= d x we have the definition of the k-th derivative of the immanant. 

3 First Expression 

We start by introducing some notation. Given a matrix A e M n (C), we will 
represent by Au\ the z-th column of A, i — {1, . . . , n}. 

Let be a natural number, 1 < k < n. Take A, X 1 , . . . ,X k e M n (C), 
and tx,...,tk k unknowns. Let Qk,n be the set of strictly increasing maps 
{1, . . . , k} — > {1, . . . , n) and Gk, n the set of increasing maps. 

We will denote by A(a; X 1 , . . . , X k ) the matrix of order n obtained from 
A replacing the a(j) column of A by the a(j) column of XK 



D> m{ XK...,X'): =1 —^- i 
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Theorem 3.1 (First expression). For every 1 < k < n, 

D k d x (A)(X\...,X k )=J2 E d x A(a;X^\...,X^). 

In particular, 

D k d x {A)(X,...,X) = k\ d x A{a;X,...,X). 

Proof. Just like in the case of the first derivative, we have that D k d x (A) (X 1 , . . . , X k ) 
is the coefficient of t\ . . . in the expansion of the polynomial d x {A + tiX l + 
. . . + tkX k ). Using the linearity of the immanant function in each column, 
we obtain the desired equality. □ 

Proof. One simply has to observe that each summand in A X (A; X 1 , . . . , X k ) 
appears (n — k)\ times: once we fix a permutation of the matrices X 1 , . . . , X k , 
these summands correspond to the possible permutations of the n—k matrices 
equal to A. □ 

As an immediate consequence of this result, we can obtain another for- 
mula for the derivative of order k of the immanant map. This generalizes 
formula (26) in [7]. 

Proposition 3.2 (First expression, rewritten). 

D k d x (A)(X\...,X h ) = j-^A x (A-,X\...,X k ). (3) 

4 Laplace Expansion for Immanants 

We start by generalizing the Laplace expansion, known for the determinant 
and the permanent, to all immanants. This expansion was proved first for 
the determinant and the same arguments used can be used to prove the 
correspondent expansion for the permanent. These can be found in [TT] and 
in [H. 

The similarity of proofs is due to the fact that the determinant and the 
permanent of a direct sum of matrices is just the product of the determinants, 
or the permanents, of the summand blocks. However, if x is an y other 
irreducible character, there is no clear general relation between the immanant 
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of A and the immanant of any submatrix of A. So the Laplace expansion 
formula for any immanant is a little more complicated. 

Let a G Q k , n , and 

S ka = {a E S n : <t({1, ...,k}) = Ima}. 
Fact 1. Ifa^/3 then S ka n S fc/J = 0. 

Suppose a G S ka fl S^. Then <t({1, . . . , k}) = Ima = Im/3, and thus 
Ima = Im/3. Since a,j3 G <5fc,n, h follows that a — j3. 

Fact 2. 

= [J Sfc^ 

If 7r G 5 n then 7r({l, ...,k}) = {j u . . .,j k }. Suppose that j\ < ■ ■ ■ < jk- 
Let 7 G S n such that 7(2) = jj. Therefore 7r G Sfc 7 and 5 n C S ka - 

The other inclusion is trivial. 

n' 

Fact 3. \{S ka : a G Q fc ,„}| = \Q k , n \ 



k\(n-k)V 



Fact 4. IfcrE S ka then cr({k + 1, . . . , n}) = Im a ( the complement of Im a ). 

Suppose that / G {k + 1, . . . , n} and a (I) = ji G Ima. We have a G S ka 
so we can find % G {ji, . . . ,j k } such that a(i) = ji. But % ^ I. This is a 
contradiction, because a G S^, and so cr is injective. 

Using the previous facts, we can conclude that for every a G S k/3 , the 
value 

n 

Yl x(v)Y[at*(t) 

v£S k/3 t=i 

does not depend on the values of the following entries of the matrix A: 



I. The first k rows of A and the n — k columns of A in Im (3. 
II. Rows k + 1, . . . , n of A and the /c columns of A with index in Im j3. 
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We now denote by A{id|/3} = (a*-) the matrix of order n obtained by 
replacing every entry in I and II by zeros. We then have 

n 

E xW)X\at«(t)=d x {A{i&\l5})- 

ves kfj t=i 

On the other hand, 

n 

d x (A{id\(3}) = J2^U a U) 
aes„ t=i 

n 

CTe U aeQfein s ka t=i 

n 

Now take a ^ (5 and er e Sfc a . Then a* CT ( t ) = 0, because at least one of 

t=i 

the factors is zero. So 

n 

yi ^ n a Ut) = °. 

^e5 fcct t=i 

for every a G \ {/?}■ 

Moreover, for A{id|/5}, if a E S kfj then a* a ^ = a t<T ( t y So rf x (A{id|/3}) = 

n 

Ttaes^xi^Yla^t). Finally, 

n 

d x( A ) = x(v)Y[at*(t) 

o-es„ t=l 

= £ rf x(^{id|«})- (4) 

Definition 4.1. Let X be a square matrix of order n, k a natural number, 
l</c<neo;,/3G <5fe,n- We denote X[a|/3] as the square matrix of 
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order k obtained from X by picking the rows a(l), . . . , a(k) and the columns 

m. 

We denote X(a\/3) as the square matrix of order n — k obtained from X 
by deleting the rows a(l), . . . , a(k) and the columns /3(1), . . . , (3{k). 

Let a, (3 G Qk,ni an d take A a square matrix of order k and B a square 
matrix of order n — k. Denote by a be the unique element of Q n -k,n with 
Ima = Ima. 

We now define 

as a square matrix of order n such that 

• Xij = if i G Im a and j ^ Im /3; 

• = if i Im a and j G Im j3; 

• Xjj = a, Q --i(j)/3-i(j) if * G Ima and j G Im/3; 

• = if z g Ima and j G" Im/3. 

In a sense, we pick rows a and /3 of A and rows a and (3 of 5. If 
a = /3 = (1, . . . , k), this is the usual direct sum of A[l, . . . , k\l, . . . , k] and 

m\ /•• i /••). 

So, according to our definitions, it is easy to verify that 
X{id|/3} = X[id|/3]0X(id|/3). 

id 

Applying the previous equality and (@J it follows that 

d x (X)= d x (X[id\0\®X(id\l3))= d x(X{id\P}). 

P£Qk,n id|/3 /3eQ fe ,„ 

For a fixed a G Qk, n , using similar arguments, we get the Laplace expan- 
sion for immanants: 

d x (X)= d x (X[a\0\®X(a\0))= £ d x(X{a\P}). (5) 
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Example 4.2. Let A be a matrix of order 4. Let k = 2 and a = id. Then, 
we have 



d x (A) 





1 On 

^11 


<7 10 
t6 lz 


o 


\ 




1 n 11 





n 1 Q 

^13 


\ 


= d x 


Gh9 
^ ±z 









%3 




«34 




aio 





«32 


zo 






&34 




Vo 





«43 


044/ 






042 





^44/ 




/ an 








a M ^ 




/o 


«12 


ai3 




+ d x 










^33 


a24 



+ d x 




^31 


«22 



«23 






C-34 




\o 


042 


043 


oy 




\a 4 i 








0-44/ 




/o 


ai2 









(0 





ai3 


a M \ 


+ d x 





«22 





033 


024 




+ d x 




031 




032 


«23 




024 






\a 4 i 





«43 






\a 4 i 


a 42 








We list some properties of the matrix 

Proposition 4.3. Let a, (3, a', (3' £ Qk,n- Then we have 

X{a\p}[a\/3] = X[a\P] and X{a\P}(a\P) = X(a\f3). 

If P ^ P' , then both matrices X {ot\fi}[ot\fi'\ and X{a\P}(a\P') have a zero 
column. If a 7^ a', then both matrices X{a\P}[a'\P] and X{a\P}(a'\P) have 
a zero row. 

Proof. These are consequences of the definitions. □ 

We can now check that this formula generalises the known Laplace for- 
mulas for the determinant and the permanent (see [11] and [Mj). If X — £ 
then d £ = det. For a £ Qk, n , denote |a| = a(l) + . . . + a{k). Fixing a £ Qk,n 
and applying the previous properties we have: 

detX = ^ det(X{a|/3}) 

= E (-!) M det(^{«|/3}[«| T ]) det(^{«|/3}(«| T )) 

/3eQ fc , n 7GQfc,n 

= (-l)' a l Y, (-l) m det(X{a\P}[a\P}) det(X {a\P}(a\P)) 

= (-l)W (-l) l/3| det(X[a|/3])det(X(a|/3)). (6) 
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This is exactly the expression of the the Laplace expansion for determi- 
nants. With similar arguments we can prove the same result for the Laplace 
expansion of the permanent. 



5 Second Expression 

We now wish to separate the entries of the matrices X 1 , . . . ,X k from the 
entries of A. It is easy to express the determinant of a direct sum in terms of 
determinants of direct summands, and the same happens with the permanent. 
With immanants, the best one can do is use formula (J5]), which is what we 
do in this second expression. 

Take X 1 , . . . , X k complex matrices of order n, and, for a G S&, G Qk,n- 
Denoting by the zero matrix of order n, we define 

XI = 0(/3;X <T(1) ,...,X CT(fc) ), 

the matrix whose 0(p)-th column is equal to X\ph,)} an d the remaining columns 
are zero, for 1 < p < k. 

Theorem 5.1 (Second Formula). 

D k d x A(X\...,X k )=J2 E WW](BM*\13)), 

in particular 

D k d x A(X,...,X) = k\ J2 d x (X[a\/3}Q)A(a\t3)) 

Proof. We have proved that 

D k d x A(X\...,X k )=J2 E d x A(P;X^\...,X^). 

By the Laplace expansion for immanants, for every G Qk,n, we have that 
d x A{(3;X a ^,...,X a(k ^ = 
= d x (A{P;X°W,...,X°W){a\l3}) 

= d x (A(ftX« 1 \...,X«Q)[a\0\@A{p ] X« 1 \...^ 
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Now we just notice that 

A(P;X^\...,X^)[a\P]=X%[a\P] 

and 

A(P;X^\...,X^)(a\P) = A(a\P). 
This concludes the proof of the formula. □ 

6 Formulas for the k-th derivative for the ra- 
th x-symmetric tensor power 

In this section we wish to establish a formula for the k-th derivative of the 
X-symmetric tensor power of a matrix. Before we can do this, we need quite 
a bit of definitions, including the very definition of this matrix. 

We start with some classical results that can be found in [131 chapter 6]. 
Let x be an irreducible character of S m and 

where id stands for the identity element of S m . K x is a linear operator 
on ® m V, and it is also an orthoprojector. Is is called a symmetriser map. 
The range of K x is called the symmetry class of tensors associated with the 
irreducible character \ an d it is represented by V x = K x (® m V). 

It is well known that the alternating character x{°~) = £ a (sign of the 
permutation a) leads to the symmetry class A m V] on the other hand the 
principal character x{°) = 1 leads to the symmetry class V m V . 

Given a symmetriser map K x , we denote 

vi* v 2 * ■ ■ ■ * v m = K x (vi ®v 2 ®...® v m ). 

These vectors belong to V x and are called decomposable symmetrised tensors. 

Let r m>n be the set of all maps from the set {l,...,m} into the set 
{1, . . . , n}. This set can also be identified with the collection of multiindices 
{(ii, . . . ,i m ) : ij < n}. If a G r mj7l , this correspondence associates to a the 
m-tuple (a(l), . . . , a{m)). In the set T m ^ n we will consider the lexicographic 
order. The set 

{aa : a G S m } C T 

m,n 
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is the orbit of a. The group S m acts on T m n by the action (a, a) — > aa^ 1 
where a G S m and a G r m n . The stabiliser of a is the subgroup of S m 
defined as 

G a = {a G S m : aa = a}. 
Let {ei, . . . , e„} be an orthonormal basis of the vector space V. Then 

{e® = e Q( i) <g> e a(2 ) ® . . . <8> e a(m) : a G r m> „} 

is a basis of the m-th tensor power of V. So the set of all decomposable 
symmetrised tensors spans V x . However, this set need not to be a basis of 
V x , because its elements might not be linearly independent, some of them 
may even be zero. Let 

n = n x = {a g r m ,„ : ^ x {a) ± o}. (7) 

creG a 

With simple calculations, we can conclude that 

ml 

creGa 

So the nonzero decomposable symmetrised tensors are {e* : a G fi}. Now, 
let A be the system of distinct representatives for the quocient set V m>n /S m , 
constructed by choosing the first element in each orbit, for the lexicographic 
order of indices. It is easy to check that A C G mi „, where G m ^ n is the set of 
all increasing sequences of T m ^ n . Let 

A = Ann. 

It can be proved that the set {e* : a G A} is linearly independent. We have 
already seen that the set {e* : a G spans V x , so there is a set A, such 
that ACACfland 

£' := {e* : a G A}, (9) 

is a basis for V x . It is also known that this basis is orthogonal if x is a linear 
character. 

In general, if x does not have degree one, there are no known orthonormal 
bases of V x (S m ) formed by decomposable symmetrised tensors. Let £ = (v a : 
a G A) be the orthonormal basis of the m-th x-symmetric tensor power of the 
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vector space V obtained by applying the Gram-Schmidt orthonormalization 
procedure to £'. Let B be the t x t change of basis matrix, from £ to 
£' = (e* : a G A). This means that for each a G A, 

7GA 

We note that this matrix B does not depend on the choice of the or- 
thonormal basis of V, since the set A is independent of the vectors, and has 
a natural order (the lexicographic order), which the basis £ inherits. More- 
over, the Gram-Schmidt process only depends on the numbers (e*,e^), and, 
by p. 163], these are given by formula 

ueSm t=i 

Hence, they only depend on the values of (e*, tj) = 6ij and are thus indepen- 
dent of the vectors themselves. 

It is known that if T is a linear operator on V, then V x is invariant for its 
m-th fold tensor power ® m T. Thus, the x-symmetric power T, denoted by 
K X (T) is defined as the restriction of ® m T to V x . There is a close connection 
of this x-symmetric tensor power of T and the immanant, as the following 
result already shows (this is a rephrasing of [HI p. 230]). 

Theorem 6.1. Suppose x is an irreducible character of the group S m . Let 
E = {ei, . . . , e n } be an orthonormal basis of the inner product space V . Let 
T G L(V, V) be the unique linear operator such that M(T, E) = A. 
Ifa,(3e r m>n , then 

(K x (T)(e* a ),e* p ) = ^d x (A T [a\(3}). (10) 
As we came across the immanant of a transpose, we note that d x (A T ) = 
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d x (A), where x is also an irreducible character of S m , defined as x( a ) := x{°)- 

m 

d x (A T ) = £ x(*)U( AT Ui) 

creS m i=l 
m 

= Yl x ^ \\ a o(i)i (i = <r _1 0')) 



CTGSm 1=1 

m 
rrt 

= xO- -1 )II a Mj) 

res m i=i 

m 

= x(r)Y[a jT{j) 

r€S m j=l 

= d Y {A) 

X tern o mesmo A que o x? [E' so' curiosidade] 



(a = r- 1 ) 



Now we want to define K X (A), the m-th x-symmetric tensor power of the 
matrix A. A natural way to do this is to fix an orthonormal basis E in V, 
and consider the linear endomorphism T such that A = M(T, E). Then the 
basis £ is orthonormal and one can define 

K X {A) :=M(K X (T),S) 

The matrix K X {A) has order t— |A|, with \Q m , n \ < t. 

It is important to notice that this matrix does not depend on the choice 
of the orthonormal basis E of V. This is an immediate consequence of the 
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formula (IT0|) : for a, j3 6 A, the (a, /3) entry of if x (A) is 

(K x (T)vp,v a ) = ( b ~fP K x( T ) e v b s a e* s ) 

7,<5eA 

= Yl b ^ b ^( K x( T ) e v e *s) 

7,5eA 



^ E b^b Sa d x (A T m) 



mi 

7,5eA 



X(id) \- , r- 



m! 

7,<5eA 

= 6 7 ^(A[5| 7 ]). 

7,5eA 

This definition admits, as special cases, the m-th compound and the m-th 
induced power of a matrix, as defined in [13, p. 236]. The matrix Kx(A) is 
called the induced matrix in [131 P- 235], in the case when the character has 
degree one. 

Denote by imm x (A) the square matrix with rows and columns indexed 
by A, whose (7,5) entry is <i x (A[7|<5]) (one could call the elements of this 
matrix immanantal minors indexed by A, the usual minors are gotten by 
considering the alternating character, in which case A = Q m _ n ). With this 
definition, we can rewrite the previous equation as 

K x (A) = ^B*imm Y (A)B. (11) 

Finally, denote by miximm x (X 1 , . . . , X n ) the square matrix with rows 
and columns indexed by A, whose (7,5) entry is A X (X 1 [7|5], . . . , X ra [ 7 |(!>]), 
so that miximm x (yl, . . . , A) = imm x (v4). We use the same shorthand as with 
the mixed immanant: for k < n, 

miximm x (A; X 1 , . . . , X k ) := miximm x (A, . . . ,A, X 1 , . . . X k ) 

Before our main formula, we recall a general result about derivatives, 
which follows from the definition. 
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Lemma 6.2. If f and g are two maps such that f o g is well defined, with g 
linear, then 

D\f o g){A){X\ ...,X k ) = D k f{g{A)){g{X l ), . . . , g{X k )). 
Theorem 6.3. According to our previous notation, we have 

D k K x (A)(X 1 ,...,X k ) = ^^ k y B* mhtimm x (A;X 1 ,...,X k )B 

and, using the notation we have already established, the (a, (3) entry of this 
matrix is 



X(id) 



E^E E Mm-y]r\p\r]®A[6m P \r)). 



ml 

7,(5GA &£S k p,TE.Qh,m p\t 

Proof. Notice that the map A h-> A[5|7] is linear, so we can apply the previous 
lemma in order to compute the derivatives of the entries of the matrix K X (A). 
The (a, 0) entry of the fc-th derivative of the m-th ^-symmetric tensor power 
of A, i.e., the (a,/3) entry of the matrix D k K x (A)(X 1 , . . . , X k ) is: 

^ £ b^b7 a D k d x (A[5\ 1 ])(X'[5\ 1 ] : ...,X k [5\ 1 ]y 

7,<5eA 

To abbreviate notation, for fixed 7, 5 G A, we will write C := A[5|7], and 
Z l := X*[5|7], i — 1, . . . , k. Using formula (J3J), we get 

D k ddA[6\i\)(X 1 [6\i\,...,X k [6\i\) = D k d x (C)(Z\...,Z k ) 

A x (C;Z\...,Z k ). 



ml 



(m — k)\ 

So the (a, 13) entry of D k K x (A)(X\ ...,X k )is 

Eb-tpbta- '—A X (C; Z 1 ,..., Z k ) = 
(m — k)l 

7,<5eA 
7,<5eA 
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According to the definition of miximm^( J 4; X 1 , . . . , X k ), we have 



D k KJA)(X 1 ,...,X k ) = X^~L- B* miximm ¥ M; X 1 ,..., X k )B. 
•° [m-k)\ 



This establishes the first formula. For the entries of the matrix, we use 
the formula in theorem 15.11 

D%(C)(Z\...,Z k )=J2 E d x{KW]@C{p\r)l 
recalling that 

where denotes the zero matrix of order m. 

So, the (at, /3) entry of the k-th derivative of K X (A) is: 

^E^E E dxiZMr]Q)C{p\T)) = 

lib. 

7,(5eA 0"£jfc P,T£Qk,m P T 

^E^aE E <h(xmz\p\T](BA[6\i\(p\T)) 

This concludes our proof. □ 

The formula obtained for the higher order derivatives of K X (A) (X 1 , . . . , X k ) 
generalises the expressions obtained by Bhatia, Jain and Grover ([7], [8]). We 
will prove this for the derivative of the m-th compound, establishing that, 
from the formula in theorem I6.3[ one can establish formula (2.5) in [ID], 
from which the main formula for the derivative of the m-th compound of A 
is obtained. 

So we take x — £ then 

K x (A)(X 1 ,...,X k ) = A m (A)(X\...,X k ). 

In this case A = Q m ^ n and the basis {e^ : a G Q m ,n} is orthogonal and it is 
easy to see (by direct computation or using formula (|S])) that every vector 
has norm 1 / y/mJ.. So the matrix B of order is diagonal and its diagonal 
entries are equal to y/mJ.. 

We now notice two properties that we will use in our computations: 
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I. For any square matrices X G M k (C), Y G M n _ fc (C) and functions 

a,/3 e Qk, n , 

detX0y = (-l) H+l/3| detXdetY. 

a\/3 

This is a consequence of formula (jSJ). We again notice that if 7 ^ /3, the 
matrices 

(X0y)[a| 7 ] and (l©F)(a| 7 ) 

a|/8 a|/3 

have a zero column. Now, using the Laplace expansion for the determinant 
along a, 

de tX0r = (-l)H £ (_i)l7l de t((X0y)H 7 ])det((X0y)( a | T )) 

a\/3 76Qfc,n a|/3 ct|/9 

= (-1) det((X Y)[a|/3]) det((X y)(«|/3)) 

"1/3 a|j8 

= (-i)W+l0ldetXdetY. 

II. For a, /3 G Q m ,n and p, r G <3fc, m , we have 

£ det(X[a|/3]^|>|T]) = fclAp^HflbM, ■ ■ ■ ,X fc [«|/3][p|r]) 

To check this, consider the columns of the matrices involved. Remember that 

X[a\0]Z = 0(T,X°V>[a\P], . . . , X^[a\/3}). 

For given a G S k and j G [A;], we have: 
entry of X[a\/3^[p\r] = entry (p(z), r(j)) of X[a\P]Z 

= entry (p(z), r(j)) of X°®[a\P] T{j) 
= entry (z,j) of X°^[a\P\W]- 

Therefore, 

X[a|/3£[p|r] = [X^WMrh] ■ ■ ■ X^[a\(3}[p\r] [k] ) 

and the matrices that appear in the first sum are the same as the ones that 
appear in the mixed discriminant. 
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We are now ready to prove the result. If we replace in theorem 16.31 
d x = det, we have that that the (a, f3) entry of D k A m (A)(X 1 , . . . , X k ) is 

^ E 6 >^E E de WH>M0^l7](pk)) 
= ^[ m! E E det(X[a|^[p|r]0A[a|/3](p|r)) 
= E E (-l) lpl+|T| det(A[a|/3](p|r))det(X[a|/3]-[p|r]) 
= fc! E (-l) H+|T| det(A[a|/3](p|r))A(X 1 [«|/3][p|r],...,X fc [a|/3][p|r]) 

We denoted by A(Bi, B n ) the mixed discriminant. The formula we 
obtained is formula (2.5) in [TO] , if you take into account that in this paper 
the roles of the letters k and m are interchanged. 

Using similar arguments one can obtain the formula for the fc-th derivative 
of V m (A)(X 1 ,...,X k ) in [8]. 
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